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Abstract 

We consider the capacitated domination problem, which models a service-requirement as- 
signing scenario and which is also a generalization of the dominating set problem. In this 
problem, we are given a graph with three parameters defined on the vertex set, which are cost, 
capacity, and demand. The objective of this problem is to compute a demand assignment of 
least cost, such that the demand of each vertex is fully-assigned to some of its closed neighbours 
without exceeding the amount of capacity they provide. In this paper, we provide the first 
constant factor approximation for this problem on planar graphs, based on a new perspective 
on the hierarchical structure of outer-planar graphs. We believe that this new perspective and 
technique can be applied to other capacitated covering problems to help tackle vertices of large 
degrees. 

1 Introduction 

For decades, Dominating Set problem has been one of the most fundamental and well-known prob- 
lems in both graph theory and combinatorial optimization. Given a graph G = {V, E) and an 
integer k, Dominating Set asks for a subset D <^V whose cardinality does not exceed k such that 
every vertex in the graph either belongs to this set or has a neighbour which does. As this problem 
is known to be NP-hard, approximation algorithms have been proposed in the literature p!tllO|[T2]. 

A series of study on capacitated covering problem was initiated by Guha et ah, [9], which 
addressed the capacitated vertex cover problem from a scenario of Glycomolecule ID (GMID) 
placement. Several follow-up papers have appeared since then, studying both this topic and related 
variations [^l [71 [8]. These problems are also closely related to work on the capacitated facility 
location problem, which has drawn a lot of attention since 1990s. See [31 116| . 

Motivated by a general service-requirement assignment scenario, Kao et al., [131 [ll] considered 
a generalization of the dominating set problem called Capacitated Domination, which is defined 
as follows. Let G = (V, E) be a graph with three non-negative parameters defined on each vertex 
u £ V, referred to as the cost, the capacity, and the demand, further denoted by w{u), c{u), and 
d{u), respectively. The demand of a vertex stands for the amount of service it requires from its 
adjacent vertices, including the vertex itself, while the capacity of a vertex represents the amount 
of service each multiplicity (copy) of that vertex can provide. 

'This work was supported in part by the National Science Council, Taipei 10622, Taiwan, under the Grants 
NSC99-2911-I-002-055-2, NSC98-2221-E-001-007-MY3, and NSC98-2221-E-001-008-MY3. 
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By a demand assignment function / we mean a function which maps pairs of vertices to non- 
negative real numbers. Intuitively, f{u,v) denotes the amount of demand of u that is assigned to 
V. We use Ng{v) to denote the set of neighbours of a vertex v 

Definition 1 (feasible demand assignment function). A demand assignment function / is said to 
be feasible if ^^ueNdv] fi'^-''^) — d{v), for each v £ V, where A''G'[f] = Ng{v) U {v} denotes the 
neighbours of v unions v itself. 

Given a demand assignment function /, the corresponding capacitated dominating multi-set 
is defined as follows. For each vertex v £ V, the multiplicity of v in 'D(f) is defined to be 



c{v) 



Xf{v) 



. The cost of the assignment function /, denoted if (/), is defined to be 



Definition 2 (Capacitated Domination Problem). Given a graph G = {V,E) with cost, capacity, 
and demand defined on each vertex, the capacitated domination problem asks for a feasible demand 
assignment function / such that w{f) is minimized. 

For this problem, Kao et al., ^1^, presented a (A-l-l)-approximation for general graphs, where A 
is the maximum vertex degree of the graph, and a polynomial time approximation scheme for trees, 
which they proved to be NP-hard. In a following work [13], they provided more approximation 
algorithms and complexity results for this problem. On the other hand, Dom et al., [6] considered a 
variation of this problem where the number of multiplicities available at each vertex is limited and 
proved the W[l ]-h.ardness when parameterized by treewidth and solution size. Cygan et al., [5], 
made an attempt toward the exact solution and presented an 0(1.89"") algorithm when each vertex 
has unit demand. This result was further improved by Liedloff et al., jl5] . 



Our Contributions We provide the first constant factor approximation algorithms for the ca- 
pacitated domination problem on planar graphs. This result can be considered a break-through 
with respect to the pseudo-polynomial time approximations given in [13], which is based on a dy- 
namic programming on graphs of bounded treewidth. The approach used in [13] stems from the fact 
that vertices of large degrees will fail most of the techniques that transform a pseudo-polynomial 
time dynamic programming algorithm into approximations, i.e., the error accumulated at vertices 
of large degrees could not be bounded. 

In this work, we tackle this problem using a new approach. Specifically, we give a new perspec- 
tive toward the hierarchical structure of outer-planar graphs, which enables us to further tackle 
vertices of large degrees. Then we analyse both the primal and the dual linear programs of this 
problem to obtain the claimed result. We believe that the approach we provided in this paper can 
be applied to other capacitated covering problems to help tackle vertices of large degrees as well. 



2 Preliminary 

We assume that all the graphs considered in this paper are simple and undirected. Let G = {V, E) 
be a graph. We denote the number of vertices, \V\, by n. The set of neighbors of a vertex v G V 
is denoted by Ng{v) = {u : {u,v) G E}. The closed neighborhood of v G ^ is denoted by 
-^gM = ^g{v) U {v}. We use degciv) and degdv] to denote the cardinality of Ng{v) and Ng[v], 
respectively. The subscript G in Ng[v] and degciv] will be omitted when there is no confusion. 

A planar embedding of a graph G is a drawing of G in the plane such that the edges intersect 
only at their endpoints. A graph is said to be planar if it has a planar embedding. An outer-planar 
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subject to 






} , J[u,v) - d[u) > U, 


u £ V 




v(^N[u\ 






c{u)x{u) - J2 ^ 0> 






v£N[u] 






d{v)x{u) — f{v, u) > 0, 






f{u,v) >0, x{u) eZ+U{0}, 


u,v £ V 


(1) 



graph is a graph which adopts a planar embedding such that all the vertices lie on a fixed circle, 
and all the edges are straight lines drawn inside the circle. For k > 1, /c-outerplanar graphs are 
defined as follows. A graph is 1-outerplanar if and only if it is outer-planar. For k > 1, a graph 
is called A:-outerplanar if it has a planar embedding such that the removal of the vertices on the 
unbounded face results in a (fc — l)-outerplanar graph. 

An integer linear program (ILP) for capacitated domination is given in The first in- 

equality ensures the feasibility of the demand assignment function / required in Definition [T] In 
the second inequality, we model the multiplicity function x as defined. The third constraint, 
d[v)x{u) — f{v, u) > 0, which seems unnecessary in the problem formulation, is required to bound 
the integrality gap between the optimal solution of this ILP and that of its relaxation. To see 
that this additional constraint does not alter the optimality of any optimal solution, we have the 
following lemma. 

Lemma 1. Let / be an arbitrary optimal demand assignment function. We have d{v) ■ Xf{u) — 
f{v, u) >0 for alluGV and v £ N[u]. 

Proof of Lemma^ Without loss of generality, we may assume that d{v) > f{v,u). For otherwise, 
we set f{v,u) to be d{v) and the resulting assignment would be feasible and the cost can only be 
better. If Xf{u) = 0, then we have f{v,u) = by definition, and this inequality holds trivially. 
Otherwise, if Xf{u) > 1, then d{v) ■ Xf{u) — f{v, u) > f{v, u) ■ {xf{u) — 1) > 0. □ 

However, without this constraint, the integrality gap can be arbitrarily large. This is illustrated 
by the following example. Let a > 1 be an arbitrary constant, and Tia) be an n-vertex star, where 
each vertex has unit demand and unit cost. The capacity of the central vertex is set to be n, which 
is sufficient to cover the demand of the entire graph, while the capacity of each of remaining n — 1 
petal vertices is set to be an. 

Lemma 2. Without the additional constraint d{v)x(u) — f{v,u) > 0, the integrality gap of the 
ILP on r(a) is a, where a > 1 is an arbitrary constant. 

Proof of Lemma The optimal dominating set consists of a single multiplicity of the central vertex 
with unit cost, while the optimal fractional solution is formed by spending ^ multiplicity at a 
petal vertex for each unit demand from the vertices of this graph, making an overall cost of ^ and 
therefore an arbitrarily large integrality gap. □ 

Indeed, with the additional constraint applied, we can refrain from unreasonably assigning a 
small amount of demand to any vertex in any fractional solution. Take a petal vertex, say v, from 
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T(a) as example, given that d{v) = 1 and f{v,v) = 1, this constraint would force x{v) to be at 
least 1, which prevents the aforementioned situation from being optimal. 

For the rest of this paper, for any graph G, we denote the optimal values to the integer linear 
program ^ and to its relaxation by OPT{G) and OPTf{G), respectively. Note that OPTf{G) < 
OPT{G). 



3 Constant Approximation for Outer-planar Graphs 



Without loss of generality, we assume that the graphs are connected. Otherwise we simply apply 
the algorithm to each of the connected component separately. In the following, we first classify the 
outer-planar graphs into a class of graphs called general-ladders and sh ow h ow the corresponding 
general-ladder representation can be extracted in O(nlog^n) time in q3.1 Then we consider in 



^3.2 and ^3.3 both the primal and the dual programs of the relaxation of ([T]) to further reduce a 
given general-ladder and obtain a constant factor approximation. We analyse the algorithm in ^3.4 



and extend our result to planar graphs in ^3.5 





Figure 1: (a) A general-ladder with anchor c. (b) A 2-outerplanar graph which fails to be a 
general-ladder, (c) The subdivision formed by a vertex u in an outer-planar embedding. 



3.1 The Structure 

First we define the notation which we will use later on. By a total order of a set we mean that each 
pair of elements in the set can be compared, and therefore an ascending order of the elements is 
well-defined. Let P = {vi,V2, . . . ,Vk) be a path. We say that P is an ordered path if a total order 
vi ^ V2 ^ ■ ■ ■ ^ Vk or Vk ^ Vk^i -<...-< wi is defined on the set of vertices. 

Definition 3 (General-Ladder). A graph G = {V, E) is said to be a general-ladder if a total order 
on the set of vertices is defined, and G is composed of a set of layers {Ci,C2-, ■ ■ ■ ,/^fc}i where each 
layer is a collection of subpaths of an ordered path such that the following holds. The top layer, £i, 
consists of a single vertex, which is referred to as the anchor, and for each 1 < j < k and u,v £ Cj, 
we have (1) N[u] C Cj-i U Cj U £j+i, and (2) u -< v implies raa,Xp(,N[u]nCj+i P =4 mingg7v[i;]n£,+i 1- 

Note that each layer in a general-ladder consists of a set of ordered paths which are possibly 
connected only to vertices in the neighbouring layers. See Fig. [T] (a). Although the definition of 
general-ladders captures the essence and simplicity of an ordered hierarchical structure, there are 
planar graphs which fall outside this framework. See also Fig. [T] (b) . 

In the following, we state and argue that every outerplanar graph meets the requirements of 
a general-ladder. We assume that an outer-planar embedding for any outer-planar graph is given 
as well. Otherwise we apply the O(nlog^n) algorithm provided by Bose [2] to compute such an 
embedding. 
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Let G = (y, E) be an outer-planar graph, u G y be an arbitrary vertex, and £ be an outer- 
planar embedding of G. We fix u to be the smallest element and define a total order on the vertices 
of G according to their orders of appearances on the outer face of if in a counter-clockwise order. 
For convenience, we label the vertices such that u = vi and vi < V2 < < ■ ■ ■ < Vn- 

Let N{u) = I^TTi , Uttji • • • ) ''^7rrfe5(„) I denote the neighbours of u such that ^ v.,^^ -< . . . ^ 
'^■^degiu) ■ -^('") divides the set of vertices except u into deg{u)+l subsets, namely. So = {^2,^3, . . . , Vt^i}, 
Si = {t;^,,w^^+i,. . . for 1 < i < deg{u), and Sdeg(u) = {^'vrd.gc^) ; ^Tr^esM+i' • • • ' ^n}- See 

Fig. [1] (c) for an illustration. 





• 




^ ^ / 




\ / 

> c 

\ _ ^ 

^ \ 
1 s 

miiibgflj 6 max„gLj a 


Figure 


2: A contradiction led by 


minfteH, 


h < maxaez,. a. 



For any < i < j < deg{u), v^^^-^ -< p ^ v^^, and f7rj_i ~< q ^ Vn^, there is no edge connecting 
p and q. Otherwise it will result in a crossing with the edge {u, v.^.), contradicting to the fact that 
is a planar embedding. 

For 1 < i < deg{u), we partition Si into two sets Li and Ri as follows. Let ds- denote the distance 
function defined on the induced subgraph of Si. Let Li = : v G Si,ds. [v-j^^^v) < ds^ {v,Vni+i)} 
and Ri = 5j\Lj. 

Lemma 3. We have max^g^. a -< min^g/j- b for all 1 < i < deg{u). 

Proof of Lemma\^ For convenience, let a = max^gj;,. a and h = min;,g/j. h. Since LiCi Ri = (p, we 
have a ^ b. Assume that b ~< a. Recall that in an outer-planar embedding, the vertices are placed 
on a circle and the edges are drawn as straight lines. Since £ is an outer-planar embedding, the 
shortest path from o to Vt^- must intersect with the shortest path from b to w-TTi+i- Let c be the 
vertex for which the two paths meet. Since Li and Ri form a partition of Si, either c £ Li 01 c £ Ri. 

If c G Li, then d^. (c, WttJ < (c, f 77^+1) by definition, which implies that ds^{b,v.,^.) < 
ds- (6, Utt-^-^) , a contradiction to the fact that b £ Ri. On the other hand, if c G Ri, then 
dsi (c, ^^ttJ > dg. (c, UjTi+i), and we have ds- (a, fjrj > ds^ (a, f^r.^-^), a contradiction to the fact 
a £ Li. In both cases, we have a contradiction. Therefore we have a ~< b. □ 

Let £{v) = dG{u,v) and ii{v) = min for any I <i < deg{u) and 

V £ Si. Observe that i{v) = ii{v) + 1, for any 1 < i < deg{u) and v £ Si. Now consider the set 
of the edges connecting Lj and Ri. Note that, this is exactly the set of edges connecting vertices 
on the shortest path between Vt^. and max^g^^ a and vertices on the shortest path between v-^^^^-^ 
and minfjgR. b. We have the following lemma, which states that, when the vertices are classified by 
their distances to u, these edges can only connect vertices between neighbouring sets and do not 
form any crossing. See also Fig. [3] 

Lemma 4. For any edge {p,q), p £ Li, q £ Ri, connecting Li and Ri, we have 
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Figure 3: Partition of Si into Lj and 



• \Kp)-Kq)\ < 1, and 

• $ edge (r, s), (r, s) ^ {p, q), r G Li, s G Ri, such that £{r) = £{q) and i{p) = £{s). 

Proof of Lemma^ The first half of the lemma follows from the definition of £. If \£{p) — £{q)\ > 1, 
without loss of generality, suppose that £{p) > £{q) + 1, by going through {p, q) then following the 
shortest path from q to u, we find a shorter path for p, which is a contradiction. The second half 
follows from the fact that f is a planar embedding. □ 

Below we present our structural lemma, which states that, when the vertices are classified by 
their distances to u, these edges can only connect vertices between neighbouring sets and do not 
form any crossing. 

Lemma 5. Any outer-planar graph G = {V, E) together with an arbitrary vertex u ^ V \s a, 
general-ladder anchored at u, where the set of vertices in each layer are classified by their distances 
to the anchor u. 

Proof of Lemma^ We prove by induction on the number of vertices of G. First, an isolated vertex 
is a single-layer general-ladder. For non-trivial graphs, let 5o,5i, . . . ,Saf,g{u) be the subsets defined 
as above. By assumption, the induced subgraphs of Sq and 5^eg(u) are general-ladders with anchors 
v-jr-^ and f7r^gg(„), respectively. Furthermore, the layers are classified by ^ — 1. That is, vertex v 
belongs to layer £{v) — 1. Similarly, the induced subgraphs of Li and Ri are also general-ladders 
with anchors and ^77^+1 whose layers are classified by £i. 

Now we argue that these general-ladders can be arranged properly to form a single general- 
ladder with anchor u and layers classified by £. Since there is no edge connecting p and q for any 
p,q with f7r,;_i -< p -< w-TTi and Vt^._-^ ^ q ^ v-j^^, < i < j < deg{u), we only need to consider the 
edges connecting vertices between Lj and Ri. By Lemma |4j when the general-ladders Lj and Ri are 
hung over Vt^. and ^tt^+i, respectively, the edges between them connect exactly only vertices from 
adjacent layers and do not form any crossing. Therefore, it constitute as a single general-ladder 
together with u and the lemma follows. □ 

Extracting the General-ladder Let G = {V, E) be the input outer-planar graph and u he 
an arbitrary vertex. We identify the corresponding general-ladder as follows. 

Compute the shortest distance of each vertex v ^ V to denoted by £{v). Let M = 
Tnax.y^y £{v) . We create M + 1 empty queues, layer{0),layer{l), . . . ,layer{M), which will be 
used to maintain the set of layers. Retrieve an outer-planar embedding of G and traverse the outer 
face, starting from u, in a counter-clockwise order. For each vertex v visited, we attach v to the 
end of layer{£{v)). 
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Theorem 6. Given an outer-planar graph Q and its outer-planar embedding, we can compute in 
linear time a general-ladder representation for G. 

Proof of Theorem^ Since the number of edges in a planar graph is linear in the number of vertices, 
the shortest-path tree computation takes linear time. The traversal of the outer face also takes 
linear time. □ 

For the rest of this paper we will denote the layers of this particular general-ladder representation 
hy Co, Ci, . . . ,Cm- The following additional structural property comes from the outer-planarity of 
Q and our construction scheme. 

Lemma 7. For any < i < M and v G Ci, we have \N{v) D Ci-i\ < 2. Moreover, if v has 
two neighbours in Ci, say, vi and V2 with vi ■< v < V2, then there is an edge joining vi (and V2, 
respectively) and each neighbouring vertex of v in Ci-i that is smaller (larger) than v. 

Proof of Lemma^ First, since the layers are classified by the distances to the anchor u, if | A^(f ) n Ci- 
3, then consider the shortest paths from vertices in N{v) n Ci-i to u. At least one vertex would be 
surrounded by other two paths, contradicting the fact that Q is an outer-planar graph. 

The second part is obtained from a similar argument. Let v' G N{v) H Ci-i be a neighbour of 
V in Ci-i. If v' is not joined to either vi or V2, then consider the shortest paths from u to vi, v' , 
and V2, respectively, v' would be a vertex in the interior, which is a contradiction. □ 



The Decomposition The idea behind this decomposition is to help reduce the dependency 
between vertices of large degrees and their neighbours such that further techniques can be applied. 
To this end, we tackle the demands of vertices from every three layers separately. 

For each < i < 3, let T^j = IJ^-^Q^sj+j. Let Qi = (Vi,£'j) consist of the induced subgraph of 
IZi and the set of edges connecting vertices in TZi to their neighbours. Formally, Vi = U^)67^i N[v\ 
and Ei = yj^^-j^. y]ui^N[v] ^iu,v). In addition, we set d{v) = for all v G Qi\TZi. Other parameters 
remain unchanged. 

Lemma 8. Let /j, < i < 3, be an optimal demand assignment function for Qi. The assignment 
function / = Ylo<i<3 fi ^ 3-approximation of Q. 

Proof of Lemma^ First, for any vertex v ^ V, the demand of v is considered in Qi for some 
< i < 3 and therefore is assigned by the assignment function fi. Since we take the union of the 
three assignments, it is a feasible assignment to the entire graph Q. 

Since the demand of each vertex in ^j, < i < 3, is no more than that of in the original 
graph G, any feasible solution to G will also serve as a feasible solution to Qi. Therefore we have 
OPT{Qj) < OPT{G), for < j < 3, and the lemma follows. □ 

3.2 Removing More Edges 

We describe an approach to further simplifying the graphs Qi, for < i < 3. Given any feasible 
demand assignment for Qi, we can properly reassign the demand of a vertex to a constant number 
of neighbours while the increase in terms of fractional cost remains bounded. 

For each v G TZi, we sort the closed neighbours of v according to their cost in ascending order 
such that w {^^{l)) < w {TTy{2)) < . . . < w {^^{deglv])), where tt^ : {1, 2, ... , (ieg'[u]} — )■ N[v] is an 
injective function. For convenience, we set -Ky{deg[v] + 1) = 0. Suppose that v G Ci. We identify 
the following four vertices. 
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• Let j„, 1 < jv < deg[v], be the smallest integer such that c (tTvUv)) > d{v). If c {tTvUv)) < d{v) 
for all 1 < j < de(7[t'], then we let jv = deg[v] + 1. 

• Let kv, 1 < kv < jv, be the integer such that and w (iTvikv)) /c{TTv{kv)) is minimized, kv is 
defined only when jv > 1. 

• Let pv = max„g7v[i,]n£^_i u and 
Qv = max„gjvHn£^+i u. 




Figure 4: Incident edges of a vertex v £ Ci to be 
kept. 

Intuitively, iTvijv) is the first vertex in the sorted list whose capacity is greater than d{v), and 
■Kv{kv) is the vertex with best cost-capacity ratio among the first jv — 1 vertices, pv and qv are the 
rightmost neighbour of v in layer C^^i and respectively. 

We will omit the function tt^, and use jv, kv to denote vr^,(j„), 7rt,(/ct,) without confusion. The 
reduced graph %i is defined as follows. Denote the set of neighbours to be disconnected from 
V by R{v) = N[v\\ {Ci U {jv U /c^ Up^ U Qv}), and let Ui = Qi\ [jvan, V}ueR{v) {(i{u,v)}. Roughly 
speaking, in graph T-Li we remove the edges which connect vertices in T^j, say u , to vertices not in 
T^j, except possibly for jv, kv, Pv, and q^. See Fig.|4j Note that, although our reassigning argument 
applies to arbitrary graphs, only when two vertices are unimportant to each other can we remove 
the edge between them. 

Lemma 9. In the subgraph Hi, we have 

• For each v ^ TZi, at most one incident edge of v which was previously in Qi will be removed. 

• For each v £ TZi, the degree of v in Tii is upper-bounded by 6. 

. OPTfiTii) <2-0PTf{gi) 

Proof of Lemma^ For the first part, let v ^ TZi be a vertex and denote S = N[v] CiTZi the set of 
neighbours of v that are in TZi. By the definition of general-ladders, for any u £ S, u ^ max^g^ a, 
we have either pu = v or qu = v, since v serves as the rightmost neighbour of u. Therefore, by our 
approach, only the edge between v and max^g^ a will possibly be removed. 

For the second part, for any v G TZi, v has at most two neighbours in the same layer, since each 
layer is a subgraph of an ordered path. We have removed all the edges connecting v to vertices not 
in Ck, except for at most four vertices, jv, kv, Pv, and qv Therefore deg{v) < 6. 

Now we prove the third part of this lemma. Let fg. be an optimal demand assignment for 
Qi, and xg^ be the corresponding multiplicity function. Note that, from the second and the third 
inequalities of ([T]), for each v £ V and u £ N[v], we have 

fg^iv,u) fgSv,u) \ , 
d{v) ' c{u) J ■ ^ ' 



xg^iu) > max 
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Maximize > a[u)yu 
















subject to 
























Vu ^ -\- (jv,Ul 


V G N[u], 






VU > 0, Zu> 0, Qy^u > 0, 


V £ N[u], 




(3) 



For each v G TZi and u G R{v) such that fgi{v,u) ^ 0, we modify this assignment as follows. If 
T^v^i"^) ^ jv, then we assign it to jv instead of to u. Otherwise, we assign it to k^. That is, depending 
on whether 7r~^(u) > j^, we raise either fgi{v,jv) or fgi{v,ky) by the amount of fgiiv,u) and then 
set fg-{v,u) to be zero. Note that, after this reassignment, the modified assignment function fg. 
will be a feasible assignment for Tii as well. 

In order to cope with this change, xg.^j^) or xg.{ky) might have to be raised as well until both 
the second and the third inequalities are valid again. If 'k~^{u) > jy, then xg^{jy) is raised by 
at most max{fg.{v,u)/d{v), fg^{v,u)/c{jy)}, which is equal to fg^{v,u)/d{v), since c(j^) > d{v). 
Hence the total cost will be raised by at most 

■ < ■ ma. I^f^. < . 

d{v) { d(v) c{u) J 

by equation Q and the fact that w{jy) < w{u). Similarly, if tt:j^{u) < jy, the cost is raised by at 
most 



d{v) ' c{ky) 

w{ky) ^ w{u) 



C{ky) 



c{u) 



since we have -^jpry < by definition of ky and equation (2). 

In both cases, the extra cost required by this specific demand reassignment between v and u is 
bounded by w{u) ■ xg^{u). Since, by Lemmajoj we have at most one such pair for each u ^ TZi, the 
overall cost is at most doubled and this lemma follows. □ 



We also remark that, although OPTf{T-Li) is bounded in terms of OPTf{Qi), an a-approximation 
for Tii is not necessarily a 2Q;-approximation for Qi. That is, having an approximation A with 
OPT{A) < a ■ OPT{rLi) does not imply that OPT{A) < 2a ■ OPT{g,), for OPTiHi) could be 
strictly larger than OPTf{T-Li). Instead, to obtain our claimed result, an approximation with a 
stronger bound, in terms of OPTf{T-li), is desired. 



3.3 Greedy Charging Scheme 

In this section, we show how we can further approximate the optimal solution for the reduced 
graph Tii by a primal-dual charging argument. We apply a technique from [Hj to obtain a feasible 
solution for the dual program of the relaxation of ([T]), which is given in ([s]) and for which we 
will further provide a sophisticated analysis on how the cost of each multiplicity we spent can be 
distributed to each unit demand it covered in a careful way. Thanks to the additional structural 
property provided in Lemma [TJ we can further tighten the approximation ratio. 
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We first describe an approacli to obtaining a feasible solution to (|3| and how a corresponding 
feasible demand assignment can be found. Note that any feasible solution to ^ will serve as a 
lower bound to any feasible solution of ([T]) by the linear program duality. 

During the process, we will maintain a vertex subset, V^, which contains the set of vertices 
with non-zero unassigned demand. For each u £ V, let d'^(ii) = J2veN[u]nV'i> ^(^) denote the 
amount of unassigned demand from the closed neighbours of u. We distinguish between two cases. 
If c(u) < d'^{u), then we say that u is heavily-loaded. Otherwise, u is lightly-loaded. During the 
process, some heavily-loaded vertices might turn into lightly-loaded due to the demand assignments 
of its closed neighbours. For each of these vertices, say v, we will maintain a vertex subset D*(v), 
which contains the set of unassigned vertices in A'^[f;]ny^ when v is about to fall into lightly-loaded. 
For other vertices, D*(v) is defined to be an empty set. 

Initially, = {u : n G Ci,d{u) / 0} and all the dual variables are set to be zero. We increase 
the dual variable i/u simultaneously, for each u G V^. To maintain the dual feasibility, as we 
increase yu, we have to raise either or gv,u, for each v E N[u]. If v is heavily-loaded, then we 
raise Zv Otherwise, we raise gv,u- Note that, during this process, for each vertex u that has a 
closed neighbour in V^, the left-hand side of the inequality c{u)zu + Ylv(^N[u] d{v)gu,v < w{u) is 
constantly raising. As soon as one of the inequalities c{u)zu + J2veN[u] ^iy)9u,v < w{u) is met with 
equality (saturated) for some vertex u G 1/, we perform the following operations. 

If u is lightly-loaded, we assign all the unassigned demand from N[u] n to u. In this case, 
there are still c{u) — d'^{u) units of capacity free at u. We assign the unassigned demand from D*{u), 
if there is any, to u until either all the demand from D*[u) is assigned or all the free capacity in 
u is used. On the other hand, if u is heavily-loaded, we mark it as heavy and delay the demand 
assignment from its closed neighbours. 

Then we set Qu = N[u]r\ and remove N[u] from . Note that, due to the definition of 
d*^, even when u is heavily-loaded, we still update d'^{p) for each p ^ V with N[p] n N[u\ ^ 0, if 
needed, as if the demand was assigned. During the above operation, some heavily-loaded vertices 
might turn into lightly- loaded due to the demand assignments (or simply due to the update of d'^). 
For each of these vertices, say v, we set D*{v) = N[v] n {V* U Q„). Intuitively, D*{v) contains the 
set of unassigned vertices from N[v] n when v is about to fall into lightly-loaded. 

This process is continued until = (j). For those vertices which are marked as heavy, we iterate 
over them according to their chronological order of being saturated and assign at this moment all 
the remaining unassigned demand from their closed neighbours to them. A high-level description 
of this algorithm is given in Fig. [5] 

Let /* : V X V ^ U {0} denote the resulting demand assignment function, and x* : V ^ 

U {0} denotes the corresponding multiplicity function. The following lemma bounds the cost of 
the solution produced by our algorithm. 

Lemma 10. For any T-Li obtained from a general-ladder Qi, we have w{f*) < 7 • OPTf{'Hi). 

Proof of Lemma \ 1 0[ We argue in the following that, for each u £V, the cost resulted by n, which 
is w{u) ■ x*{u), can be distributed to a certain portion of the demands from the vertices in N[u] 
such that each unit demand, say from vertex v £ N[u], receives a charge of at most 7 ■ y^. 

Let M € y be a vertex with x*{u) > 0. If u has been marked as heavy, then by our scheme, 
we have gu,v ~ and y^ = Zy, for all V £ N[u]. Therefore 'w{u) = c{u) ■ Zu = c{u) ■ y^, and 
for each multiplicity of u, we need c{u) units of demand from A^[m]. If x*{u) > 1, then at least 
c{u) ■ {x*{u) — 1) units of demand are assigned to u, and by distributing the cost to them, each 
unit of demand gets charged at most twice. If x*{u) = 1, then we charge the cost to any c(n) 
units of demand that are counted in d'^{u) when u is saturated. Since u is a heavily- loaded vertex, 
d'^{u) > c{u) and there will be sufficient amount of demand to charge. 
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Algorithm Greedy- Charging 

1: i — {u:ueV, diu) ^ 0}, V* < — V. 
2: < — Et,eAfH ^(^)' fo'^ ^^'^^ ueV. 
3: w'^{u) i — w{u), for each n € 
4: Let Q be a first-in-first-out queue. 
5: while / do 

6: ry < — w'^{v)/ in.in{c{v) , d'f' {v)} , for each v G V*. 

7: ii — argm.in{ry : v G V*}. [u is the next vertex to be saturated.] 

8: w'^{v) i — w'^{v) - w'^{u), for each v eV*. 

9: 

10: if d'^{u) < c{u) then 

11: Assign the demand from N[u\ fl to u. 

12: Assign c{u) — d'^{u) amount of unassigned demand from D*{u), if there is any, to u. 
13: else 

14: Attach u to the queue Q and mark u as heavy. 
15: end if 

16: Let Su< — N[u]nV'f'Ll{u}. 

17: Remove N[u] from and update the corresponding d'^{v) for v G V. 

18: For each v E V* such that d'^{v) = 0, remove v from F*. 

19: for all vertex v becomes lightly-loaded in this iteration do 

20: D*{v)i — N[v]n{V^USu). 

21: end for 

22: end while 

23: while Qy^ (f) do 

24: Extract a vertex from the head of Q, say u. 
25: Assign the unassigned demand from N[u] to u. 
26: end while 

27: 



Figure 5: The high-level pseudo-code for the primal-dual algorithm. 

On the other hand, if u is lightly-loaded, then x* (u) = 1 and we have two cases. If J2v&N{u] "^l^) — 
c{u), then u is lightly-loaded in the beginning and we have Zy = and y„ = gu,v for each v G N[u], 
which implies w(u) = J2veN[u] '^(^) ' 9^,'" ~ "llveNlu] '^(^) " y-"- '^^^ w{u) of u can be distributed 
to all the demand from its closed neighbours, each unit demand, say from vertex v G N[u], gets a 
charge of yy. 

If ^veN[u] '^('") ^ c('u), then u is heavily-loaded in the beginning and at some point turned into 
lightly-loaded. Let Uq C D*{u) be the set of vertices whose removal from makes this change. 
By our scheme, Zy is raised in the beginning and at some point when d'^{u) is about to fall under 
c(u), we fixed Zy and start raising gy^y for v G D*{u)\Uq. Note that, we have = Zy for all 
Uq G Uq., yv = Zy + gu,v for each v G D*{u)\Uo, and gu,v = for t> G N[u]\D* {u) U Uq. Let (i^^ = 

c(^^)-E^ei?*(»)Wo ^(^)- ^^^^ ^(") = c(tx)-z„ + 5]^g^[„] d{v)-gu,y = (dl^^ + EyeD*{u)\Uo ^iv)) ■ 

Zu+^yeD*{u) d{v)-gu,v = dli^-Zu+Y.v&D*{u)\Uo d{v)-{zu+gu,v) < E^GWo dvyv+T.v&D*{u)\Uo d{v)-yy. 
Again the cost w{u) of the single multiplicity can be distributed to the demand of vertices in D*{u). 

Finally, for each unit demand, say demand d from vertex ii, consider the set of vertices Vd C N[u\ 

that has charged d. First, by our assigning scheme, consists of at most one heavily-loaded vertex. 
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If d is assigned to a heavily- loaded vertex, then, by our charging scheme, we have |V^| = 1, and d is 
charged at most twice. Otherwise, if d is assigned to a lightly-loaded vertex, then, by our charging 
scheme, each vertex in charges d at most once, disregarding heavily-loaded or lightly-loaded 
vertices. 

By the above argument and Lemmajoj we have w{f*) = J2ueV w{u)x*{u) < Xluey '^(^) 'deg[u] ■ 
Uu ^ 7 • l^uev ^(^) — 0PTf{7ii), where the last inequality follows from the linear program 
duality of Q and ([s]). □ 






Figure 6: Situations when a unit demand of u is fully-charged. 



Thanks to the structural property provided in Lemma [7j given the fact that the input graph 
is outer-planar, we can modify the algorithm slightly and further improve the bound given in the 
previous lemma. To this end, we consider the situations when a unit demand from a vertex u with 
(ie5f[ii] = 7 and argue that, either it is not fully-charged by all its closed neighbours, or we can 
modify the demand assignment, without raising the cost, to make it so. 

Lemma 11. Given the fact that Tii comes from an outerplanar graph, we can modify the algorithm 
to obtain a demand assignment function /* such that < 6 • 0PTf{7ii). 



Proof of Lemma 11. Consider any unit demand, say demand d from vertex u in Cj, and let C 
N[v\ be the set of vertices that has charged d by our charging scheme. 

First, we have |iV[u]| < 7 by Lemma [9| By our charging scheme, |A^[ti]| < 7 implies \Vd\ < 7. 
In the following, we assume |A^[ti]| = 7 and argue that either we have \ Vd\ < 7, or we can modify 
the solution such that |Vd| < 7. By assumption, ku, Pu, Qu are well-defined. Let ui,U2 € N{u) CiCj 
denote the set of neighbours of u in Cj such that ui ■< u < U2- By Lemma [7| depending on the 
layer to which j„ and belong, we have the following two cases. 



Both ju and k^ belong to /Zj+i. If two of ku, say, j„ and ku, are not joined to ui and 
U2 by any edge, then at most one of ju and k^ can charge d, since u is the only vertex with possibly 
non-zero demand in their closed neighbourhoods. When the first closed neighbour of u is saturated 
and u is removed from V^, both ju and ku will be removed from V* and will not be picked in later 
iterations. Therefore, at most one of ju and ku can charge d. 

On the other hand, if two of {ju-,ku,qu]-, say ju and qu, are joined to ui and U2, respectively, 
then we argue that at most two out of {ju,u,qu} can charge d. Indeed, ui, u, and U2 are the only 
vertices with non-zero demands in the closed neighbourhoods of {ju, u, qu}- After two of {ju, u, qu} 
is saturated, ui, u, and U2 will be removed from V^. Therefore, at most two out of {ju,u,qu} can 
charge d. See also Fig. [6] (a) and (b). 

Only one of {ju,ku} belongs to >Cj+i and the other belongs to £j i. Without loss of 
generality, we assume that ju G ^^j-i and ku G ^j+i- By Lemma [7| both ju and pu are joined 
either to ui or U2 separately. Since Hi is outerplanar, we have ju < u < pu, otherwise ju will be 
contained inside the face surrounded by pu, ui, and u, which is a contradiction. See also Fig. M (c). 
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If both ku and g„ are not joined to either ui or U2, then by a similar argument we used in 
previous case, at most one of ku and Qu can charge d. Now, suppose that, one of {ku, Qu}, say, ku, is 
joined to ui by an edge. We argue that, if both ui and ku have charged d after d has been assigned 
in a feasible solution returned by our algorithm, then we can cancel the multiplicity placed on ku 
and reassign to ui the demand which was previously assigned to ku without increasing the cost 
spent on ui. 

If ui is lightly-loaded in the beginning, then the above operation can be done without extra 
cost. Otherwise, observe that, in this case, ui must have been lightly-loaded when d is assigned 
so that it can charge d later. Moreover, ui is also in the set and not yet served, for otherwise 
ku will be removed from V* and will not be picked later. In other words, at this moment when ui 
becomes lightly- loaded, we have ui G V^, meaning that it is possible to assign the demand of ui 
to itself later without extra cost. 

On the other hand, if ui or ku is the first one to charge d, consider the relation between Qu 
and M2- If there is no edge between Qu and U2, then Qu will not be picked and will not charge d. 
Otherwise, if {qu,U2) exists in Hi, then it is a symmetric situation described in the last sequel. 

In both cases, either we have \ Vd\ < 6 or we can modify the solution returned by the algorithm 
to make \V(i\ < 6 to hold. Therefore we have w{f*) < 6 • OPTf{'Hi) as claimed. □ 



3.4 Overall Analysis 

We summarize the whole algorithm and our main theorem. Given an outer-planar graph G = {V, E), 



we use the algorithm described in ^3.1 to compute a general-ladder representation of G, followed 
by applying the decomposition to obtain three subproblems, Qq, Qi, and Q2- For each Qi, we use 
the approach described in |3.2|to further remove more edges and obtain the reduced subgraph T-Li, 



for which we apply the algorithm described in ^3.3 to obtain an approximation, which is a demand 



assignment function /j for Hi. The overall approximation, e.g., the demand assignment function 
/, for G is defined as / = Eo<j<3 /«• 

Theorem 12. Given an outerplanar graph G as an instance of capacitated domination, we can 
compute a constant factor approximation for G in 0{v?) time. 



Proof of Theorem 12 First, we argue that the procedures we describe can be done in 0{n ) time. 
It takes O(nlog'^n) time to compute an outer-planar representation of G f2]. By Theorem |6| 
computing a general-ladder representation takes linear time. The construction of Qi takes time 
linear in the number of edges, which is linear in n since G is planar. 

In the construction of the reduced graphs Hi, for each vertex v, although we use a sorted list of 
the closed neighbourhood of v to define , k^, Pv, and Qy , the sorted lists are not necessary and Tii 
can be constructed in 0{n) time by a careful implementation. This is done in a two-passes traversal 
on the set of edges of Gi as follows. In the first pass, we iterate over the set of edges to locate jy, 
Py, and Qy, for each vertex v G V. Specifically, we keep a current candidate for each vertex and 
for each edge (u, v) £ E iterated, we make an update on u and v if necessary. In the second pass, 
based on the jy computed for each v G V, we iterate over the set of edges again to locate ky. The 
whole process takes time linear in the number of edges, which is 0{n). 

In the following, we explain how the primal-dual algorithm, i.e., the algorithm presented in 
Figure [sj can be implemented to compute a feasible solution in O(n^) time. First, we traverse the 
set of edges in linear time to compute the value d^{v) for each vertex. In each iteration, the next 
vertex to be saturated, which is the one with minimum w'^{v)/ min {c{v), d'^{v)^, can be found in 
linear time. The update of w'^{v) for each v G V* described in line 8 can be done in linear time. 
When a vertex v G N[u] with non-zero demand is removed from V^, we have to update the value 



13 



d'^{v') for all v' G N[v]. By Lemma [oj the closed degree of such vertices is bounded by 7. This 
update can be done in 0(1) time. The construction of Su can be done in linear time. Since d'^{v) 
can only decrease, each vertex can turn into lightly-loaded at most once. Therefore the process time 
for these vertices is bounded in linear time. The outer-loop iterates at most 0{n) times. Therefore 
the whole algorithm runs in O(n^) time. 

The feasibility of the demand assignment function / is guaranteed by Lemma [8] and the fact 
that Tii is a subgraph of Gi. Since Hi C Q^^ the demand assignment we obtained for Tii is also a 
feasible demand assignment for Gi. Therefore, / is feasible for G. 

By the definition of /, Lemma 11 Lemma |9j and Lemma [8j we have 



(/)< ^(^*) H OPTf{n,)< 12- OPTfiG,) 

0<j<3 0<i<3 0<j<3 

<12- ^ OPT{Gi) <2,Q-0PT{G). 

0<i<3 

□ 



3.5 Extension to Planar Graphs 

We describe how our outer-planar result can be extended to obtain a constant factor approximation 
for planar graphs under a general framework due to [1]. This is done as follows. Given a planar 
graph G, we generate a planar embedding and retrieve the vertices of each level using the linear-time 
algorithm of Hopcroft and Tarjan |llj . 



Pv 




qv 



Figure 7: (a) 3-outerplanar graph, (b) Local connections with respect to a 
vertex v. Bold edges represent links in the ladder extracted from Li. Thin 
edges represent links between Li and the other two levels, Lq and L2. 



Let m be the number of levels of this embedding. Let OPT be the cost of the optimal demand 
assignment of G, and OPTj be the cost contributed by vertices at level j. For convenience, in the 
following, for j < or j > m, we refer the vertices in level j to an empty set and the corresponding 
cost OPTj is defined to be zero. 

For < i < 3, we define Gi as ^o<j<f {OPT^.j^i + OPT^.j-i+2)- Since Ylo<i<z Ci <2- OPT, 
there exist an iq with < io < 2 such that Ci^ < | • OPT. For each < j < ^ , define the graph Gj 
to be the graph induced by vertices between level 3 ■ j — iq and level 3 ■ j — zq + 2. The parameters 
of the vertices in Gj are set as follows. For those vertices who are from level 3 ■ j — iq and level 
3 ■ j — io + 2, their demands are set to be zero. The rest parameters are remained unchanged. 
Clearly, Gj is a 3-outerplanar graph, and we have X]o<j<— OPT{Gj) < OPT + | • OPT, where 
OPT{Gj) is the cost of the optimal demand assignment of Gj. 

In the following, we sketch how our algorithm for outerplanar graphs can be modified slightly 
and applied to Gj for each < j < y to obtain a constant approximation for Gj. For convenience, 
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we denote the set of vertices from the three levels of Gj by Lq, Li, and L2, respectively. See 
Fig. [7] (a). 

• Obtaining the General Ladders. For each level Lj, < i < 3, which constitutes an 
outerplanr graph by itself, we define a total order over it according to the order of appearances 
of the vertices in counter-clockwise order. The general-ladder is extracted from Li as we did 
before. Furthermore, for each vertex in the ladder, its incident edges to vertices in Lq and L2 
are also included. 

• Removing Redundant Edges. In addition to the four vertices we identified for each vertex 
V with non-zero demand, we identify two more vertices, which literally corresponds to the 
rightmost neighbours of v in level Lq and L2, respectively. See also Fig. [7](b). Then, the first 
part of Lemma [9] still holds, the degree bound provided in the second part is increased by 2, 
and the third part holds automatically without alteration. 

The rest parts of our algorithm remain unchanged. From the above argument, we have the 
following theorem. 

Theorem 13. Given a planar graph G as an instance of capacitated domination, we can compute 
a constant factor approximation for the G in polynomial time. 

4 Conclusion 

The results we provide seem to have room for further improvements. One reason is that, due to 
the flexibility of the ways the demand can be assigned, it seems not promising to come up with an 
approximation threshold. However, when the demand cannot be split, it is not difficult to prove 
a constant approximation threshold. Therefore, it would be very interesting to investigate the 
problem complexity on planar graphs. 

Second, as we have shown in §3.1[ the concept of general-ladders does not extend directly to 
fe-outerplanar graphs for k > 2. It would be interesting to formalize and extend this concept to 
/c-outerplanar graphs, for it seems helpful not only to our problem, but also to most capacitated 
covering problems as well. 
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